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Abstract 

The  results  of  mathematical  end  experimental  studies  presented  in  preceding 
reportectearty  show  toat  toe  model  A/i  accurately  describes  the  properties  of  a 
frozen  fringe  when  the  steady  "fnwto  of  an  ice  layer  oocure.  In  this  work  the 
steady  growth  of  ice-rich  frozen  s  Is  studied  by  using  Mj.  Deriving  a  traveling 
wave  sot  ution  to  the  problem,  we  tu»  found  that  the  condition  of  steady  growth 
of  Ice-rich  frozen  soil  is  uniquely  determined  by  a  set  of  two  physical  variables, 
such  as  ao  and  ctj  used  earlier,  under  given  hydraulic  conditions  and  over¬ 
burden  pressures  arta  that  the  traveling  wave  solution  converges  to  the  solu^ 
to  the  problem  of  a  steadily  growing  ice  layer  when  the  velocity  of  toe  0°C 
isotherm  relative  to  the  unfrozen  part  of  toe  soil  vanishes. 
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NOMENCLATURE 

Oq  function  defined  by  eq  63g 

U\  function  defined  by  eq63h 

b  defined  by  eq  63c 
b0  function  defined  by  eq  63d 
bj  function  defined  by  eq  63e 
b2  function  defined  by  eq  69b 
b3  function  defined  by  eq  69c 

B,  ith  constituent  of  die  mixture.  Subscripts  i  = 

1, 2,  and  3  are  used  to  denote  unfrozen  water, 
ice  and  soil  minerals,  respectively 

c  heat  capacity  of  the  mixture  defined  by  eq 
10c 

Co  defined  by  eq  41c 
Cj  heat  capacity  of  the  ith  constituent 
d  unit  of  time,  day 

dt  density  of  the  ith  constituent 

fi  mass  flux  of  the  ith  constituent  relative  to 
that  of  soil  minerals  where  i  =  1, 2 

/y  mass  flux  of  the  ith  constituent  relative  to 
that  of  soil  minerals  in  Rj  where  i  ■  1, 2  and 
j-0,1,2 

hi  heat  content  of  the  ith  constituent 

I0  function  defined  by  eq  56a 

I]  function  defined  by  eq  56b 

I2  function  defined  by  eq  59b 

k  thermal  conductivity  of  the  mixture 
Jfcj  thermal  conductivity  in  R{  where  i  =  0, 1, 2 
kn  limiting  value  of  k3  defined  by  eq  40h 
k2i  limiting  value  of  defined  by  eq  48 

Ko  hydraulic  conductivity  in  die  unfrozen  part 
of  die  soil 

Kj  empirical  function  defined  by  eq  1  where  i  = 

1,2 

Kn  limiting  value  of  Kj  as  §  approaches  nt  while 
4  is  in  Rj,  i  *  1,2 

limiting  value  of  Kj  as  §  approaches  no  while 
4  is  in  JRX,  i  -  1,2 

L  latent  heat  of  fusion  of  water,  334  J  g-1 
m  boundary  where  the  content  of  unfrozen 
water  is  negligible 

Mi  name  of  a  model  defined  in  Part  I  where  i  = 
1,2,3 


n  boundary  in  Ro  and  also  used  a  generic  mov¬ 
ing  surfa« 

ri  velocity  of  n  =  dn/dt 

n j  boundary  with  i  =  0, 1  where  «o  denotes  die 

boundary  where  T  =  0°C  and  «|  the  interface 
between  R2  and  a  frozen  fringe 

p0  gravity  term,  0.098  kPa/cm 
Pj  pressure  of  die  ith  constituent  where  i  =  1,2 
Pio  value  of  Px  at  «o 
Pln  value  of  Px  at  n 
Pn  value  of  P2  at ft 

q  heat  flux  in  the  mixture  by  conduction  defined 

byeq8b 

q*  limiting  value  o(qa&$  approaches  n2  while  4 
isinRj 

cf  limiting  value  of  q  as  $  approaches  n3  while  5 
isinR2 

ft  heat  flux  in  the  rdi  constituent  by  conduction 

r  rate  of  frost  heave 

Rq  unfrozen  part  of  the  soil 

R}  frozen  fringe 

R2  frozen  part  of  the  soil 

Rm  region  in  the  diagram  of  temperature  gradi¬ 
ents  where  an  ice  layer  melts 

Rg  region  in  die  diagram  of  temperature  gradi¬ 
ents  where  the  steady  growth  of  an  ice  layer 
occurs 

R«  boundary  between  R,  and  Ru 

Ru  region  in  the  diagram  of  temperature  gradi¬ 
ents  where  die  steady  growth  of  an  ice  layer 
does  not  occur 
s  defined  by  eq  32b 

s+  defined  by  eq  31b 

S  defined  by  eq  71e 

Sj  defined  by  eq  70d-70g  where  i  =  1,2,  -,  4 

t  time 

T  temperature  of  die  mixture 

T/  temperature  at  ft  where  i  =  0, 1 

T0  temperature  at  ft)  and  used  also  as  a  reference 
temperature 

ft  velocity  of  die  ith  constituent  where « *  1, 2, 3 

ftj  ft  in  Rj  where  i  =  l,2, 3  and;  *0,1,2 


iv 


Y  defined  by  eq  20 

Vj  V  in  Kj  where;  =  0,1,2 
lu-,  defined  by  eq  35b  where  i  =  0, 1, 2 
®ij  defined  by  eq  35a  where  i  =  0, 1  and  j  =  0, 1, 
2 

x  spatial  coordinate 
X  defined  by  eqAl 

Y  defined  by  eq  B4 
Yj  defined  by  eq  B6 
z  defined  by  eq  10b 
Zj  defined  by  eq  44c 

absolute  value  ofthe  temperature  gradient  at 
«o 

at  absolute  value  of  the  limiting  temperature 
gradient  as  ^  approaches  rq  while  J;  is  in 
defined  by  eq  47 

Po  defined  by  eq  41b 

Pi  defined  by  eq  46 

Y  constant,  1.12  MPa  “C-1 

5  thickness  of  a  frozen  fringe 

Sq  defined  by  eq  54c 
e  defined  by  eq  58b 
q  defined  by  eq  40f 
6j  volumetric  content  of  file  ith  constituent 

X|  rate  of  supply  of  mass  of  the  ith  constituent 

per  unit  volume  of  the  mixture 
rate  of  surface  supply  of  mass  of  the  ith  con¬ 
stituent  per  unit  surface  of  the  mixture  de¬ 
fined  by  eq  14a 

A  function  defined  by  eq  26f 


A\  function  defined  by  eq  23f 
p  defined  by  eq24e 
Mo  defined  by  eq  34a 
v  defined  by  eq22h 

v  defined  by  eq  B8 
V}  defined  by  eq  23i 

vi  v  at  T  =  Tj 

§  coordinate  defined  by  eq  16 
\  frost  heave  ratio  defined  by  eq  72 
*0  defined  by  eq  B2 
X}  defined  by  eq  B3 
pi  bulk  density  of  the  ith  constituent 
Pi)  Pi»nXj 

o  effective  stress  defined  by  eq  60b 

ot  defined  by  eq  59a 

$o  empirical  function  of  T  defined  by  eq  55a 

$oi  value  of  $q  at  T  =  Tj 

$i  empirical  function  of  T  defined  by  eq  55b 

$u  value  of  $i  at  T  =  Tj 

$2  empirical  function  of  T  defined  by  eq  55c 

$21  value  of  $2  at  T  =  Ti 

V  some  function  of  x  and  t 

lyl  jump  of  \*f  defined  by  eq  lib 

y*  defined  by  eq  12b 

\|T  defined  by  eq  12a 

*  superscript  used  to  indicate  the  value  of  any 

variable  evaluated  when  a  point  (oq,  Oq)  in 
the  diagram  of  temperatime  gradients  is  on 

r: 


Traveling  Wave  Solutions  to  the  Problem 
of  Quasi-Steady  Freezing  of  Soils 

YOSHISUKE  NAKANO 


INTRODUCTION 

We  will  consider  the  one-directional  steady 
growth  of  an  ice  layer.  Let  die  freezing  process 
advance  from  the  top  down  and  the  coordinate  x 
be  positive  upward  with  its  origin  fixed  at  some 
point  in  the  unfrozen  part  of  the  soil.  A  freezing 
soil  in  this  problem  may  be  considered  to  consist  of 
three  parts:  the  unfrozen  part  R&  the  frozen  fringe 
Rx  and  die  ice  layer  R2  as  shown  in  Figure  1.  The 
physical  properties  of  parts  Rq  and  R2  are  well 
understood,  but  our  knowledge  of  die  physical 
properties  and  the  dynamic  behavior  of  part  Rx 
does  not  appear  sufficient  for  engineering  applica¬ 
tions. 


T*T_ 


T.T, 


T-To-0(°C) 


Figure  1.  Steady  growth  of  ice-rich  frozen  soil. 


The  results  of  our  mathematical  and  experi¬ 
mental  study  on  die  steady  growth  condition  of  an 
ice  layer  were  presented  in  die  three  previous 
reports  (Nakano  1990,  Takeda  and  Nakano  1990, 
Nakano  and  Takeda  1991).  These  results  clearly 
show  that  die  model  Mi  accurately  describes  die 
properties  of  a  frozen  fringe  during  the  steady 
growth  of  an  ice  layer  under  negligible  overbur¬ 
den  pressure.  The  model  Mx  is  the  frozen  fringe 
where  ice  may  exist  but  does  not  grow  during  the 
steady  growth  of  an  ice  layer  and  the  mass  flux  of 
water f\  is  given  as 

(1) 

drr  dx 

where  x  is  the  space  coordinate,  and  and  K2  are 
the  properties  of  a  given  soil  that  generally  depend 
cm  temperature  T  and  die  composition  of  die  soil. 

Nakano  (1990)  has  shown  that  the  velocity 
ri0(=  driQ/dt)  of  die  frost  front  is  nonpositive  and 
vanishes  when  die  steady  growth  of  an  ice  layer 
occurs .  In  this  work  we  will  study  a  case  in  which 
no  is  a  given  negative  constant  and  the  steady 
growth  of  ice-rich  frozen  soil,  instead  of  an  ice 
layer,  takes  place.  In  such  a  case  the  frozen  fringe 
Rj  also  moves  downward  with  a  constant  speed. 
Ice  may  or  may  not  exist  in  Rv  However,  if  a 
certain  steady  distribution  of  ice  is  present  in  Rt, 
then  the  growth  of  ice  must  occur  in  R\  because 
unlike  the  case  of  a  steadily  growing  ice  layer,  die 
velocity  of  soil  partidesinR}  relative  to  nodoesnot 
vanish  when  no  <  0.  Because  of  this  we  must 
modify  die  definition  ofMt  so  thatke  may  grow  in 
Ri  when  no  <  0* 


The  objective  of  this  work  is  to  show  that  there 
exists  a  traveling  wave  solution  to  the  problem  of 
steadily  growing  ice-rich  frozen  soil  and  that  this 
solution  is  reduced  to  the  solution  to  the  problem 
of  a  steadily  growing  ice  layer  obtained  in  the 
previous  report  (Nakano  1990)  when  the  velocity 
no  vanishes.  We  will  also  show  that  the  condition 
of  a  the  steady  growth  of  ice-rich  frozen  soils 
under  given  hydraulic  conditions  and  applied 
pressures  is  uniquely  determined  by  a  set  of  two 
physical  variables,  such  as  Oq  and  ax,  used  in  the 
previous  reports. 


given  as  (Nakano  1986) 

— (ftfci)  =  ~  —  (piM)-  —  <Ji,  i=lA3  (6) 
dt  dx  dx 

where  p  x  h  { is  the  heat  content  of  the  ith  constituent 
per  unit  bulk  volume  and  qx  the  heat  flux  by  con¬ 
duction.  We  assume  that  the  constituents  are  local¬ 
ly  in  thermal  equilibrium  with  each  other,  Le.,  that 
the  constituents  have  locally  a  common  tempera¬ 
ture  T  (°C).  Under  such  an  assumption,  the  heat 
content  ht  is  given  as 


BASIC  EQUATIONS 

We  will  treat  the  soil  as  a  mixture  of  water  in  the 
liquid  phase  Bv  ice  B2  and  soil  minerals  fl3  with 
bulk  densities  pi,  P2  and  P3,  respectively.  If  dj  is  the 
density  of  the  ith  constituent,  then  the  volumetric 
content  0j  of  the  ith  constituent  is  given  as 

0i  =  Pi/4-  (2) 


h1  =  c1(T-T0)  (7a) 

h2  =  -L  +  c2(T  -  T0)  (7b) 

h3  =  c3(T-T0).  (7c) 

We  sum  up  eq  6  for  i  =  1, 2, 3  to  obtain  the  heat 
balance  equation  of  die  mixture  given  as 


It  is  clear  that  the  stun  of  0j  should  be  unity,  namely 

0x  +  02  +  03  =  1.  (3) 

We  will  assume  that  the  density  of  each  constitu¬ 
ent  of  the  mixture  remains  constant.  Thus,  the  re¬ 
sults  of  this  study  are  accurate  if  the  deformation 
of  each  constituent  is  negligibly  small,  regardless 
of  overburden  pressure.  The  dry  density  of  the  un¬ 
frozen  part  of  the  soil  is  assumed  constant  during 
the  freezing  process. 

We  will  assume  that  the  unfrozen  part  of  the  soil 
is  kept  saturated  with  water  at  all  times  by  using 
an  appropriate  device  of  water  supply.  The  bal¬ 
ance  of  mass  for  the  ith  constituent  is  given 
(Nakano  1986)  as 


(8a) 

dt  i  dx  i  dx 

where  q  is  defined  as 

(8b) 

It  is  known  that  qx  depends  on  the  bulk  density  pi# 
the  thermal  properties  of  ith  constituent,  the  tem¬ 
perature  gradient  and  the  way  in  which  the  ith 
constituent  is  distributed  in  the  mixture.  For  the 
sake  of  simplicity  we  will  approximate  q  as 

?  =  -*- -  (8c) 

dx 


|-(PW)  +  *b  *  =  1/2,3  (4) 

at  ax 

where  vx  is  the  velocity  of  the  ith  constituent  and  A* 
is  the  time  rate  of  supply  of  mass  of  the  ith  constit¬ 
uent  per  unit  volume  of  the  mixture.  It  should  be 
mentioned  that  the  summation  convention  on  in¬ 
dex  i  is  not  in  force  here,  so  that  (pjdj)  represents 
only  one  term.  Since  none  of  the  constituent  is 
involved  in  chemical  reaction,  we  have 

X.j  +  —  0  and  Xj  =  0.  (5) 


where  fcisthe  thermal  conductivity  of  the  mixture 
that  generally  depends  on  the  thermal  properties 
of  each  constituent  and  the  composition  of  the 
mixture. 

Using  eq  4,  we  reduce  eq  8a  to 

£  ft  1  ki  =  -  £  -  X  Plf,  f  *i  -  f  (9) 

i  df  i  i  dx  dx 

We  will  assume  that  q  and  L  in  eq  7a,  7b  and  7c 
do  not  depend  on  T.  Choosing  T0  to  be  0  (°C)  and 
using  eq  7a,  7b  and  7c,  we  reduce  eq  9  to: 

j-q  =  I^X2  +  z) 
dx 


The  balance  of  heat  for  the  ith  constituent  is 


(10a) 


where  z  is  defined  as 

Lz  =  -c^+(c i-C5)TX2-£psmCi^I  (10b) 
at  i  dt 

C  =  C1Pi+C2p2  +  C3p3.  (10c) 

We  will  now  consider  a  moving  surface  whose 
location  is  given  as 

x  =  n(f)  <,  m  ( t ).  (11a) 

In  a  neighborhood  of  n(t)  we  choose  two  moving 
surfaces  n~(t)  and  n+(t)  with  n~(t)  >  n(t)  >  n*(t).  The 
jump  of  a  quantity  y(x,t)  at  n(t)  is  defined  as 

M  =  V“-V+  (lib) 

where 


conservation  law  of  either  heat  or  mass  must  be 
violated  if  one  of  these  conditions  does  not  hold 
true  at  n. 

Quasi-steady  problem 

We  will  consider  a  special  case  in  which  a  frost 
front  x = «0(/)  moves  with  a  constant  velocity  mq.  In 
such  a  case  we  will  seek  a  quasi-steady  solution  to 
the  problem  described  by  eq  4  and  10a  in  the  form 
of  a  traveling  wave.  We  will  introduce  a  new 
independent  variable  %  defined  as 

\  =  x-nof  (16) 

where  «o  =  drio(t)/dt.  Using  eq  16  we  reduce  eq  4 
to 

ft  (p«  -  »o)  =  Xj ,  i  =  1,  2,  3  (17) 

ac, 


\|/  =  lim  v  (12a) 

V+=  limy  (12b) 

It  is  clear  that  lyl  =  0  if  y  is  continuous  at  n(t). 

Jump  conditions  at  n(t)  under  the  assumption 
of  a  continuous  T  are  given  (Nakano  1986)  as 

|pi®il  =  |pil  n+K(  i  =  1,  2,  3  (13a) 

M  ■  |Xpifri|  n  — |X»ipifri|  (13b) 

where  Xj  is  the  surface  supply  of  the  mass  of  the  ith 
constituent  defined  as 

Xi=  lim  I  hdx  (14a) 

«+,  rr-*rt  Jn* 

and  n  is  assumed  to  be  continuous  and  is  defined 
as 

n  =  i&  (14b) 

dt 

From  eq  5  we  obtain 

Xi  +  X2  =  0  and  X3  =  0.  (15) 

The  jump  conditions,  eq  13a  and  13b,  are  necessary 
and  sufficent  conditions  for  the  conservation  law 
of  heat  and  mass  to  hold  at  n.  In  other  words,  the 


For  the  sake  of  convenience  we  will  define  new 
dependent  variables  fx  and  f2  as 

fi  =  Pi  (»i  -  W3)  (18a) 

fl  -  P2  (»2  -  »3)  (18b) 

It  is  easy  to  see  that  (i  =  1, 2)  is  the  mass  flux  of 
either  Bx  or  B2  relative  to  the  mass  flux  of  soil 
particles.  Using  eq  18a  and  18b,  we  reduce  eq  17  to 

(Pi V)'=-/{-X2  (19a) 

(pz  V)'=-/2'+X2  (19b) 

(psV)'=0  (19c) 

where  primes  denote  differentiation  with  respect 
to  ^  and  V  is  defined  as 

V=V3  -«„.  (20) 

Using  eq  16, 18a  and  18b,  we  will  reduce  eq  10a 
and  10b  to 

<r=-(*n'=  ^  +  (21a) 

Lz  =  -(ci/i+c2f2+cV)T,  +  (ci-c2)X2T.  (21b) 

Traveling  wave  solution 

We  will  derive  a  traveling  wave  solution  that 
satisfies  the  jump  conditions,  eq  13a  and  13b,  at  nv 
and  the  balance  equations  of  mass  and  heat,  eq  19a 
19b  19c  and  21a.  We  will  assume  that  the  bound¬ 
aries  rt,«o,ni  and  m  in  Figure  1  move  with  the  same 


3 


constant  velocity,  namely: 

«  =  «o  =  »i  =  rh. 

(22a) 

The  pressure  Pt  of  water  is  assumed  to  be  a  given 
constant  at  n  as 

Pi(n)  =  Pu. 

(22b) 

According  to  Mx  (Nakano  1990),  the  mass  flux  of 
water  fn  in  Rj  is  given  as 

fn  =  -Kl^l-K2~  for  4  in  Ri 

34  34 

(22c) 

K2  /Ki  -*  y  as  fn-*0 

(22d) 

lim  Pi(4)  *P2(ni)  =  P2l 

(22e) 

tin  Hi 


where  P2  is  the  pressure  of  the  frozen  part  of  the 
soil  R2  and  y  is  a  constant  with  the  value  of  1.12 
MPa  °C~\ 

We  will  also  assume  that  the  composition  of  the 
soil  is  continuous  and  X2  vanishes  at  «o  as 

|Pil  =  0,  A-2  =  0  at  no-  (22f) 


V?  and  that  of  Vj  as  4  approaches  n0witi\^>n0by 
VT  as  shown  in  Figure  2.  Our  immediate  task  is  to 
reduce  the  number  of  unknown  variables  appear¬ 
ing  in  this  figure  by  using  the  jump  conditi<Mts  arul 
the  balance  equations  of  mass  and  heat. 

Integrating  eq  19a,  19b  and  19c  from  4  =  «f 
to  m  *,  we  obtain 

fii  +  p6V2+= fn  +  (hiV2~-  At  (23a) 

^V2+=teVT+ Ai  (23b) 

(23c) 

where 

Vi  =  v&  -  no  =  r  -  no  (23d) 


At=  I  M4  (23f) 

Jnf 

where  it  is  assumed  that  the  temperature  Tm  at  4 = 
m  is  low  so  that/^  and  pi2  vanish.  This  implies  that 
is  equal  to  the  rate  of  frost  heave  r .  Eliminating 
A?,  we  will  reduce  eq  23a,  23b  and  23c  to 


The  assumption  of  eq  22f  implies  that  the  velocity 
0;and  the  flux  of  heat  q  are  continuous  at  tiq.  We 
will  assume  that  the  movement  of  ice  relative  to 
soil  particles  is  negligibly  small  everywhere,  that 
is 

/2  =  0  in  Rj  (i  =  0,1,2).  (22g) 

As  we  discussed  (Nakano  1990),  when  file  steady 
growth  of  an  ice  layer  occurs,  the  pressure  Pj  is 
continuous  but  the  first  derivative  dPx/dx  of  Px 
may  be  discontinuous  at  n0.  Therefore,  the  bound¬ 
ary  no  is  generally  a  free  boundary  where  the  first 
derivative  dPx/d£,  may  be  discontinuous.  Finally 
we  will  assume  that  px  is  given  as 

Pi  =  P3v(T)  (22h) 


va- 

0  * 
"l2  ' 

p* . 
22 

P4‘ 

*12 

r2 

v2‘- 

P12  1 

p»* 

P-  . 

f  12 

V’ 

piV 

p2V 

P*  ’ 
K31 

*lV  A+ 

ri 

vi. 

P11  • 

P21  • 

P31  • 

fn  •  A 

V 

piV 

P«’ 

P31  ‘ 

fTt 

where  v(T)  isa  given  empirically  determined  func¬ 
tion  of  T  that  is  assumed  to  be  approximated  by  the 
equilibrium  unfrozen  water  content  at  T. 

We  will  now  denote  the  values  of  vx  and  V,  for 
instance,  in  the  part  Rfj  =  0,1,2)  by  and  Vy 
respectively.  We  will  denote  the  limiting  value  of 
Vx,  for  instance,  as  4  approaches  nx  with  4  <  «i  by 


R0  v0  ’  P10  ’  P30  ■  *10 

n  - 

Figure  2.  Variables  in  Rq,  Rj  and  R2. 
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P&V?  ■  (viPsi  +  P2j)MT +/12 

(23g) 

P&V?  -  ffaV2~ 

(23h) 

where  vt  is  defined  as: 

Vi  =  vOi). 

(23i) 

Using  eq  3,  we  obtain 

Pa  =  <#2(1  -d$  f^) 

(23j) 

(fa  =  di[i  -  (vidj1  +  <f 3 pa ] . 

(23k) 

From  eq  23f,  23i  and  23j,  we  obtain 

V}=  [l  +  vj  {d2  -df1)  fe]  V2‘  +  d2 

V12.  (23 1) 

It  is  easy  to  see  that  all  variables  at  m+  are  deter¬ 
mined  if  all  variables  at  nf  are  known.  Using  eq 
13a  at  Rj,  we  obtain 

fn  +  ptiVf =fa  +  pt2V2  +  X2 

(24a) 

pfiV?"  =  P22V2  —^>2 

(24b) 

pJiVi+=  tfavi 

(24c) 

where 

K  -  oft  -  rio . 

(24d) 

For  the  sake  of  simplicity,  we  will  introduce  a  new 
variable  p  defined  as 

Pk  =  V  P31  • 

(24e) 

Using  eq  22h  and  24e,  we  will  reduce  eq  24a,  24b 
and  24c  to 

fii  +  vip$iVt  =  fa  +  vi\ip$iVi  +  X2  (25a) 

pliV^  =  P22V2  — 

(25b) 

VT=  pa¬ 

(25c) 

using  eq  3,  we  obtain 

vi  ppnd  f1  +  (fad?  +  PPsida1  =  1 

(25d) 

Vi  p$\di  +  ptid2  +  p&d3  =  1 

(25e) 

From  eq  13b,  we  obtain 


W-XjfL+fa-cjjT,].  <250 

Integrating  eq  19a  19b  and  19c  from  Roto  £,  we 

obtain  the  following  equations  given  as 

fn  +  vp3iVi  =/io  +  pioVi”-  A  (26a) 

P21^  =  A  (26b) 

p3iVi=  psoVf  (26c) 

where 

Vi  =  Vai  -  no  (26d) 

V\  =  Vo  =  -  ri0  (26e) 

A(0=  Ml  ^n0.  (260 

.  no 

Using  eq  3,  we  obtain 

vp3i  d  i 1  +  p2i  df1  +  p3i  d 31  =  1  (26g) 

pio  df1  +  P30  df1  =  1.  (26h) 


Taking  limits  of  eq  26a,  26b  and  26c  as  \  approach¬ 
es  tip  we  obtain 

fn  +  vi  V?+  A+  =  fw  +  pio  Vo  (27a) 

p?,VT-A+  =0  (27b) 

p3iVt  =  p»V0.  (27c) 

Comparing  eq  25a,  25b  and  25c  with  eq  27a,  27b 
and  27c,  respectively,  we  obtain 

fa  +  vippsV2-=/io+pioV()- A+-A,2  (28a) 

piV2'=  A+  +  X2  (28b) 

pVZ  =  (p3o/p^)V0.  (28c) 

Let  us  assume  for  the  time  being  that  rio ,  v2  (or  T2) 
and  \2  are  given.  It  is  easy  to  see  that  the  left-hand 
terms  of  eq  28a,  28b  and  28c  contain  four  unknown 
variables  at  tt~  in  Rp  Vi,  {fa,  p,  and  fn,  while  the 
right-hand  terms  of  these  three  equations  contain 
unknown  variables  in  the  combined  region  of  Rq + 
Rv  Since  {fa  and  p  are  related  by  eq  25d,  all  the 
unknown  variables  at  nf  listed  in  Figure  2  are 
uniquely  determined  if  all  the  variables  in  Rq + R\ 
listed  in  Figure  2  are  known. 
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From  eq  27a,  27b  and  27c  we  find  that  these 
three  equations  contain  five  unknown  limiting 
values,  Vf,  p$1#  pji,/?i  and  A+.  Since  pft  and  Pm 
are  related  by  eq  25e,  we  have  actually  four  un¬ 
known  limiting  values  and  three  equations.  There¬ 
fore,  if  one  of  these  four  unknowns  is  given,  then 
all  six  unknown  limiting  values  at  n  ?  listed  in 
Figure  2  are  uniquely  determined.  Choosing  V j  to 
be  an  independent  variable,  we  will  write  all  other 
limiting  values  as  follows.  First  from  eq  27c,  we 
obtain: 

P3J  =  p3o(Vo/Vt)  (29a) 

Pi+i  =  vip3o(Vo/V?).  (29b) 

From  eq  25e,  we  obtain: 

pft  =  d2-d2 (vvii'+d Jl)p»( Vo /Vl)  .  (29c) 

From  eq  27a  and  27b,  we  obtain 

fix  -fio  +  pioVo  -  (vip^i  +  p2i)V^  (30a) 

A+  =  pftVf.  (30b) 

Substituting  in  eq  30a  by  eq  29c,  we  will  reduce 

eq30ato 

fti  =fio  +  s*  Vo  -  d2(Vt  -  Vo)  (31a) 

where  s+  is  defined  as 

s+  =  (l  -  d^d^)  (pio  -  V1P30)  ■  (31b) 

Unknown  variables,  p31,  pn,  P21//11  and  A  are 
givenby  eq  29a,  29b,  29c,  31a  and  30b,  respectively, 
in  which  superscripts  +  are  deleted  and  v2  is  re¬ 
placed  by  v.  For  instance, /u  is  given  as 


fn  -fio  +  sVQ  -  d2(Vi  -  V0)  (32a) 

where  s  is  given  as 

s  =  (l  -  (pio  -  vp3o) .  (32b) 

Using  eq  28a  and  28c,  we  will  write  unknown 
variables  at  nj  as 

V2'  =  di%  +  Vt  (33a) 

P'V^  +  dz'1^)*1  (33b) 

P32  -  MPli  (33c) 


Pi2  =  v«W>3i  (33d) 

Pb  ■  d i[l  -  (virfr1  +  dj1)  mpm]  (33e) 

(33f) 

Using  eq  31a,  we  will  reduce  eq  33f  to 

fn  =/io  +  s"V0- d2 (Vi*-  Vo)  - 12 .  (33g) 

In  actual  experiments,  p10  and  £30  are  given  as 
initial  conditions.  If  Vq,  Vj  (or  T j),  X.2,  fw  and  Vj  are 
given,  then  all  other  variables  are  uniquely  deter¬ 
mined  .  Since  p£i  is  difficult  to  measure  experimen¬ 
tally,  it  is  convenient  to  introduce  a  new  variable 
Mo  defined  as 

Mo=pftpM=(pft/pfc)(pft'P»)P-  (34a) 

Using  eq  23h,  25c  and  27c,  we  will  reduce  eq  34a  to 

Mo=  Vo/Vf .  (34b) 

Using  eq  24e  and  25c,  we  will  reduce  eq  23/  to 

Vf=  M'*[l  +  viUj"1  -  dj1)  upli]  Vt+dffo  (34c) 

Substituting  fn  in  eq  34c  by  eq  33g  and  using  eq 
33b,  we  will  reduce  eq  34c  to 

V?  =  dzVio  +  V0[l  +  (d?  -  dr1)  pio] .  (34d) 

Combining  eq  34b  and  34d,  we  obtain 

Mo  =  Vo^Vio+Vo  [l  +  (d21-di1)]  pio)  V  (34e) 

From  eq  23d  and  34d  we  obtain 

r  =  d2-1/io  +  (d2-1- df1)  pioVo.  (34f) 

Using  eq  34f,  we  will  reduce  eq  34e  to 

Mo=  Votr  +  VoT1.  (34g) 

We  will  introduce  new  variables,  and  w-y  de¬ 
fined  as 

=  i  =  1,2 and j  =  0, 1, 2  (35a) 

u>j  ~  wij  +  W2j  j  =  0, 1, 2.  (35b) 

Itis  clear  thatn>jj  is  thecontentofiheith  constituent 
in  Rj  and  that  w)  is  the  content  of  ice  and  unfrozen 
water  in  Rj.  We  will  refer  to  iej  as  the  total  water 


content.  Using  vty  we  will  reduce  eq  32a  to 


fn  =/io  +  Pxfwo  -  «>i  Wo  •  (36a) 

Using  eq  26a,  26b  and  36a,  we  obtain 

A  =  P3o(a>i-v)VV  (36b) 

We  will  now  examine  the  behavior  of/n.  From  eq 
32a  we  obtain 

fn  =/io  +  sV0  -  d2v 3i .  (37a) 

Using  eq  26c,  we  will  write  u31  as: 

V31  =  (pso  pai  - 1)  Vo .  (37b) 


Since  T  decreases  as  £  increases  from  n0  to  nx,  s  is 
positive  in  Rx  and  increases  with Jj.  It  is  anticipated 
that  p31  may  decrease  with  £  but  does  not  increase 
with  Since  Vo  2  0,  from  eq  37b  we  find  that  u31 2 
0  and  v31 2: 0  in  Rx .  Therefore,  from  eq37a  we  obtain 

fw  +  s+Vo  2  fn  2  fw  -  (38a) 

where 

»3+i  =(pw-  P31)  Vo  /p3i<  (38b) 

Ice-rich  frozen  soil 

We  will  focus  the  remainder  of  our  analysis  on 
a  special  case  in  which  the  frozen  part  of  die  soil 
contains  a  significant  amount  of  ice.  For  such  a 
case  the  mobility  of  water  in  R2  is  anticipated  to  be 
much  less  than  that  in  Rx  and  we  may  neglect  fa. 
It  follows  from  eq  34a,  34g  and  36b  that  the  values 
of  po,  V0  and  A  renuuu  small. 

Theexact  composition  of  Rx  isnotknown.  How¬ 
ever,  it  is  a  generally  accepted  view  that  p31  does 
not  change  significantly  from  px.  The  results  of 
our  analysis  on  the  data  of  Tomakomai  silt  (Na- 
kano  and  Takeda  1991)  appear  to  support  such  a 
viewpoint.  Assuming  the  existence  of  a  certain 
rule  for  p3J,  we  will  explore  probable  rules  below. 
Suppose  that  such  a  rule  is  known,  then  two  of  five 
independent  variables,  Vi  and  X&  are  uniquely 
determined  by  eq  26c  and  33g,  respectively,  when 
three  remaining  independent  variables,  V^Tj  and 
fio,  are  given  in  this  case.  Let  us  consider  first  a 
special  rule  that  p31  is  kept  constant  at  pM.  In  such 
a  case  v31  vanishes  and  eq  37a  is  reduced  to 

fn  -  fio  +  sVq  •  (39a) 


When  V0>  0,  we  obtain 

fn  >  fio  ondfn  >0,  if  Vo  >  0  and  £  >  no.  (39b) 

It  follows  from  eq  39b  that/n  is  greater  than  f1Q  and 
increases  with  This  special  case  does  not  appear 
to  be  probable  because  die  mobility  of  water  should 
not  increase  with  increasing 


Next  we  will  consider  a  special  rule  that  the 
total  water  content  aq  is  kept  constant  at  w0.  From 
eq  36a  and  33f  we  obtain: 

> 

ii 

o 

II 

(40a) 

From  eq  37a  we  obtain 

V31  =  dj1  sVo. 

(40b) 

From  eq  29a  we  obtain 

psi  =  pso(l  +  dj1  • 

(40c) 

We  will  reduce  eq  36b  to 

A  =  p3o(w'o-v)V0. 

(40d) 

It  follows  from  eq  40b  and  40c  that  u31  increases 
with  £  while  p31  decreases  with  \  for  a  given  V„. 
This  second  special  case  appears  more  probable 
than  the  first  case  because  the  mass  flux  of  water 
should  not  increase  with  increasing  We  will 
study  the  second  case  below.  The  empirical  func¬ 
tion  v(T)  in  eq  40d  is  known  to  be  an  increasing 
function  of  T  with  v(0)  =  w0.  We  will  assume  that 
v(T)  possesses  a  continuous  first  derivative. 

The  thermal  conductivity  Jq  of  R1  depends  on 
the  composition  of  Rj.  Our  experimental  data  indi¬ 
cate  that  Iq  is  a  nondecreasing  function  of  £.  We 
will  approximate  Aq  by  a  linear  function  of  %  as 

*i(5)  =  *bD+n(5-»o)]'  «i>^2«o  (40e) 

n  =  (fcn-ltb)/(5fco)20  (40f) 

8  =  *!-«(>  (40g) 

lim  ki(4)  =  kn  £  *21  (40h) 

tinRi 

where  k2\  is  die  limiting  value  of  k2  when  \  ap¬ 
proaches  ti\  while  %  is  in  R2.  Under  assumptions 
described  above  we  will  study  thermal  and  hy¬ 
draulic  fields  below. 
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Temperature  Ttlj) 

We  will  seek  solutions  T(0  to  the  balance  equa¬ 
tion  of  heat  given  by  eq  21a  in  Rq  and  Rt.  We  will 
begin  with  Rq.  Since  fx  =  /10,  /2  =  0,  V=  V0and  = 
0,  from  eq  21a  and  21b  we  obtain 


r-ft,r  =  o,  SinRo 

(41a) 

Po  =  (ci/io  +  c0V0)/ko 

(41b) 

where  ko  is  he  thermal  conductivity  of  Rq  and  Cq  is 
defined  as 

Co  =  clPl0  +  c3p30  • 

(41c) 

Since  T(no)  =  0°C,  integrating  eq  41a,  we  obtain 

7t0=aopb1{l  -exp[-po(no  -01) 

(42a) 

rfe)=-cxoexp[-po(no  -0] 

(42b) 

where  Oq  is  defined  as 

ao  =  -!r(«o). 

(42c) 

Next  we  will  seek  a  solution  in  Rx.  First  we  will 
rewrite  eq  21a  and  21b  by  using  eq  26a,  26b  and 
26c.  Since  f\  =  /10,/2  =  0  and  A^  =  A'in  this  case,  we 
will  reduce  eq  21b  to 

Lz  =  -(ci/io  +  cV1)r  +  (c1-c2)AT.  (43a) 

Using  eq  26a,  26b  and  26c,  we  obtain 

cVx  =  -  (cj -Cj)A  +  c0V0 .  (43b) 

Using  eq  43b,  we  will  reduce  eq  43a  to 

Lz  =  -(Ci/io  +  cbV0)r  +  (cl-c2)(A7)'.  (43c) 

Using  eq  43c,  we  will  reduce  eq  21a  to 

q'  =  -*^r  +  (Cl-C2)(AT)'  +  LA'.  (44a) 

Integrating  eq  44a  from  %  =  ixq  to  Ij,  we  obtain 

q  =  -k{T  =  kfptq- zxT  +  LA  (44b) 

Zi  =  *oPo-(ci-C2)A  (44c) 

where  kx  is  die  thermal  conductivity  of  Rx  and  A  is 
given  as 

A  =  (it>0-v)p3oV0-  (44d) 


Using  eq  41b  and  44d,  we  will  reduce  eq  44c  to 

zi  =  Ci/10  +  P30V0  lc3  +  c2w o  +  (c,  -  cj)v].  (44e) 

Neglecting  the  last  term  in  eq  44e,  we  will  reduce 
eq44bto 

-^r^JfcoOo-MiT+LA.  (45) 

where  is  defined  as 

Pi  =  XoVl/lO  +  P3Q  V'o(C3+C2»o)].  (46) 

We  will  reduce  the  jump  condition  (eq  25f)  to  a 
somewhat  more  convenient  form  by  using  eq  44b. 
We  will  write  the  limiting  value  q~  when  \  ap¬ 
proaches  «i  while  £  is  in  R2  as 


q-  =  *2icti 

(47) 

foi  =  lim  1:2(0 

(48) 

$-»«i 

Sinlfe 

where  Jt2  is  the  thermal  conductivity  of  R2anda1  is 
the  absolute  value  of  the  limiting  temperature 
gradient  as  £  approaches  n2  while  £  is  in  R2.  Using 
eq  44b  and  47,  we  will  reduce  eq  25f  to 

*2iai  ~  +  ^bPo^i  =  {fio  +  A+)  [L  +  (ct  -  Cj)Tx J 

(49) 

A+  =  (w0  -  Vi)p3oVo .  (50) 

As  shown  in  Appendix  A,  eq  45  has  a  unique  and 
decreasing  solution  for  nx  £  £2  wq.  For  a  special  case 
in  which  the  following  condition  holds  true, 

t\5<1  and  ^S<1.  (51) 

We  found  hat  he  condition,  eq  51,  holds  true 
when  he  steady  growth  of  an  ice  layer  occurs 
(Takeda  andNakano  1990).  When  eq51  holds  true, 
we  may  use  an  approximation  (Nakano  1990)  giv¬ 
en  as 

=  1  +  no).  (52) 

We  obtain  an  approximate  solution  (App.B)  given 
as 

7t5)  =  KoteiV -(m  -  itoit^v  pi)(§  -  no )  (53a) 
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r(^X  »  -«» + *o(v  +  x[l  Piv) 

-  ftUi-1**!1  frvjfc  -  no).  (53b) 

It  is  easy  to  see  from  eq  45  that  the  temperature  T(t,) 
is  uniquely  determined  if  V^/KyOo  and  ^  are  given. 
Hence,  the  temperature  Tj  at  nx  is  determined  by 
Viv/io,  (Xq  and  5.  Suppose  that  Vq,  fl0>  Oq  and  8  are 
given,  that  Tj  is  determined  by  eq  45.  Once  Tj  is 
known,  then  cq  is  determined  by  eq  49.  As  we 
described  in  the  preceding  section,  all  variables 
listed  in  Figure  2  are  determined  by  Vq,  Tj  and  /10 
if  p31  obeys  a  certain  known  rule  in  Rj.  This  implies 
that  four  independent  variables  must  be  given  in 
order  to  specify  the  condition  of  freezing.  We  will 
choose  (Xq,  (Xj,/j0  and  V0  to  be  independent  vari¬ 
ables. 

Pressure  Pj(£) 

When  the  mass  flux  of  water  is  given  by  eq  1 ,  we 
have  found  (Nakano 1990)  that  the  following  equa¬ 
tions  hold  true  in  Rx: 

Pio  =  Pi„  ~  [(fio  /Xo)  +  po]5o  (54a) 

|'«l  rlt  1 

K?K2TdZ,  -/J  K?d%  (54b) 

no  Jno 

where  P10  =  Pi  (no),  Pi„  =  Pi  (n), 
n  =  some  point  in  Rq 
Ro  =  hydraulic  conductivity  in  Rq 
p0  =  gravity  term  that  is  equal  to  the  den¬ 
sity  dj  multiplied  by  foe  gravitation¬ 
al  acceleration. 

8q  is  defined  as 


T,*0 


IT1  (Kxo/KMki/ko)dT  T<  0 


(55b) 


*j(7)  = 


Ti  =  0 

Tri  (v /woUio  /Ri)  (fcx  /ko)dT  T  <  0 
0  (55c) 


i 


where  Kxo  and  Xgj are  the  limiting  values  of  Xj  and 
X2,  respectively,  as  %  approaches  tiq  while  £  in  Rj. 

Choosing  T  as  an  independent  variable,  we  will 
write  foe  two  integrations  in  eq  54b  as 


Jo* 


hs 


xr1  X2  Td%  =  X,'o  X20  Ti4>oi 


(56a) 

(56b) 


where  fyjj  =  ^(Tj).  We  will  write  eq  54b  as 


p2i  =  pio  “ /o -/iq/i- 


(57) 


Using  eq  55b  and  55c,  we  will  reduce  eq  56b  to 

(App. n) 

h  -  -(itiXio)  1[<|>iiTi(l  -e)+a>oWoiti1^2iT,i] 

(58a) 


£  =  Pi  k!1 


-T,+ 


(58b) 


8q  =  n0  -  n  >  0. 


(54c) 


We  will  assume  that  P21,  Pla  and  8q  are  given. 

In  order  to  reduce  eq  54b  to  a  simpler  form,  we 
will  introduce  foe  following  three  dimensionless 
quantities 


4»o(T)  > 


Ti  =  0 


1 


Tj  I  (Xio /Ki) (X2 /X20)  dT  T<  0 


(55a) 


where  <>„  =  <t>i(T2)  and  fcj  =  (Tj) . 

Using  eq  56a  and  58a,  we  will  write  eq  57  as 

°i  -  p2i  “  pio  =  ~  (59a) 

h  -  Xjo  X20  <poi  -  (« 1  Xio ) 1 

/»[♦«  (1  -  e)  +  u>0*onfl92i  ] .  (59b) 

Using  eq  54a,  we  obtain 

<Ti  =  or + p08o + SoKq1  fio  (60a) 

o  =  P21  -  pu  (60b) 
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-  Ti  =  (o  +  po8o  +  SoXo'/io)/  h  (60c) 

Since  the  composition  of  the  freezing  soil  is  as¬ 
sumed  to  be  continuous  at  %  we  may  expect  that 
the  limiting  value  K10  and  Xq  should  be  equal 

Xo  =  X10.  (61a) 

Neglecting  the  gravity  term,  we  will  reduce  eq  60c 
to 

-  Ti  =  (or  +  80K0V10)/  h-  (61b) 

When  fXQ  vanishes,  from  eq  61b  we  obtain 

a  =  -(Km  /Xo)  <|»oi  T\,  if /io  =  0.  (61c) 

Hie  generalized  Clausius-Clapeyron  equation 
(Edlefsen  and  Anderson  1943),  which  was  proven 
empirically  by  Radd  and  Oertle  (1973),  is  given  as 

o  =  —  yTj ,  if  yjo  —  0.  (62a) 

Comparing  eq  61c  with  eq  62a,  we  obtain 

r^o/Woj,  if/io  =  0.  (62b) 

It  follows  from  eq  22d  that  eq  62b  holds  true  and 
that  we  have 


Y=X2o/Xo 

(62c) 

<fcn  =  l. 

if/10  =  0. 

(62d) 

It  should  be  noted  that  eq  62c  should  hold  true 
regardless  of  fw.  Using  eq  62c,  we  reduce  eq  59b  to 

h  =  Yfon  -  («iXo)_1 

ho  foil  (1  -  e)  +  loonoJti^i].  (62e) 

For  a  special  case  in  which  o  is  negligibly  small,  eq 
61b  is  reduced  to 

-Ti  =  (5oKo1/io)//2.  (62f) 

At  the  end  of  the  preceding  section  we  had  four 
independent  variables,  Oq,  ctj,/i0  and  V0.  Since  we 
have  derived  another  equation,  eq  60c,  we  now 
have  three  independent  variables,  Oq,  oq  and  V0. 
We  will  derive  one  more  equation  below  in  order 
to  reduce  the  number  of  independent  variables  to 
two. 


Condition  of  steady  growth 

Since  the  mass  flux  is  given  by  eq  22c,  the  flux/10 
in  a  neighborhood  of  »q  is  given  as 

ho  =  -  Xu  P[(n?)  -  X21 T '(«,+)  (63a) 

whereXn,  X21  and  T\nf)  are  the  limit¬ 
ing  values  of  Xi,  X2  and  T\\)  ,  respective¬ 

ly,  as  £  approaches  tq  while  f;  is  in  Rj.  From  eq  53b 
we  obtain 


r(n,+)  =  -  aob 

(63b) 

b  -  (l  +  qSjf1  (60  +  M) 

(63c) 

b=a  o1  Ui-*ovi-*oK^piVi) 

(63d) 

b\  =  a*  ft  (ki  -  icouj1  frv,). 

(63e) 

Similarly  from  eq  53a  we  obt 

Ti  =  -ao(®o  +  «i8) 

(63f) 

aos-a^Ttoiij’v, 

(63g) 

fl^a^Oii-xon^viPi) 

(63h) 

where  vj  =  v(Ti) . 

Using  eq  63b,  we  will  reduce  eq  63u  to 

/10  =  -  Xu  P{  (nt )  +  MC21  ao . 

(64a) 

Neglecting  terms  representing  sensitive  heat,  we 
will  reduce  eq  49  to 

Mo +  ¥10  =  *21«1  -  PxWv’o  -  Vl). 

(64b) 

Now  we  will  recall  a  special  case  studied  (Nakano 
1990)  where  V0  vanishes  and  the  steady  growth  of 
an  ice  layer  occurs.  In  such  a  case.  Kg  vanishes  and 
eq  63c  is  reduced  to 

6  =  (l  +  q5)-1(l  +  p16) 

(65a) 

Pi  =  -  fro  Wio . 

(65b) 

Also  eq  64b  is  reduced  to 

fro<Xo  +  Lfio  =  fr2i«i  • 

(65c) 

It  was  found  (Nakano 1990)  that  the  steady  growth 
of  an  ice  layer  occurs  under  die  conditions  given  as 

(*21  /kofai  >  Oq  (66a) 
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Figure  3.  Temperature  gradients  dj  and  Oq. 

P{  (ni  )  >  0.  (66b) 

Using  eq  64a  and  65c  and  combining  eq  66a  and 
66b,  we  obtain 

(k2i/k0)a1  >  do  >  k2i(ko  +  IM21)_1oq .  (66c) 

The  region  Rs  in  Figure  3  satisfies  eq  66c  and  the 
steady  growth  of  an  ice  layer  occurs  in  R,.  Line  R»*  * 
in  Figure  3  is  given  as 

<*o  =  (^2i/^b)®i*  (67a) 

Line  R»**  is  the  boundary  between  R*  and  R„, 
where  an  ice  layer  melts.  The  boundary  R!  in 
Figure  3  is  given  as 

a0  =  *21  ( ko  +  Lb* Kn )  ai  (67b) 

where  superscripts  *  are  used  to  indicate  the  value 
ofany  variable  whenapoint(a1,a0)  belongs  to  R? . 
Since  b*  and  Kh  generally  depend  on  (Xq  and  oq, 
the  boundary  R?  between  R,  and  R^,  where  the 
steady  growth  of  an  ice  layer  does  not  occur,  is  a 
curve  stemming  from  the  origin.  From  eq  66b  we 
obtain 


P{  (n{ )  >  0 

in  R. 

(67c) 

P[{nt)  =  0 

on  Ri. 

(67d) 

Now  we  will  examine  the  case  in  which  V<j  ® 
positive  and  the  steady  growth  of  ice-rich  frozen 
soil  occurs.  First  we  will  show  that  the  necessary 
condition  for  fire  steady  growth  is  given  as 


Pt'(irf<)*0.  (68) 

Suppose  that  eq  68  does  not  hold  true.  Since 
Pi'  (I)  <  0  in  a  neighborhood  of  nj  in  R},  there 
exists  a  point  in  this  neighborhood  such  that 
Pi<y>p2i-  Also  P2i  >  Pm  because  Pzi  20and/io 
>  0.  Since  P](£)  is  continuous  in  Rv  there  exists  at 
least  one  point  $2  with  >  ^2  >  no  such  that  P1(^2) 
=  P2i  .  This  implies  that  another  ice-rich  frozen  soil 
layer  can  exist  in  Rj.  This  is  obviously  contradicto¬ 
ry  to  our  assumption  that  there  is  no  ice-rich  frozen 
soil  layer  in  Ry 

When  eq  68  holds  true,  from  eq  64a  and  64b  we 
obtain 


ao  2  *2i  f»2 1  Ol  -  *3 

(69a) 

where  b2  and  b3  are  defined  as 

1^2  ~  kg  +  LbK21 

(69b) 

b3  -  bfpxVoUwo  -  Vl'). 

(69c) 

In  Figure  3  we  will  draw  curve  1  given  as 

oto  =  *21  b2  ai  -  63. 

(69d) 

Since  V0  >  0,  curve  1  must  be  in  R^  and  converges 
to  R»*  when  V0  approaches  zero.  When  VQ  vanishes 
and  the  steady  growth  of  an  ice  layer  occurs,  a  line 
of  constant  fw  is  parallel  to  Ri*,  such  as  broken 
line  2  in  Figure  3.  It  follows  from  eq  64b  that  line  2 
is  still  fire  line  of  constant  f10  except  in  R„  where 
line  2  is  the  line  of  constant  fl0  +  P30V0 (“to  ~  vi)> 

It  follows  from  eq  69d  that  the  distance  between 
curve  1  and  Rf  increases  with  increasing  V0.  From 
eq  34g  we  find  that  the  ice  content  in  R2  decreases 
with  increasing  V0.  The  condition  eq  69a  implies 
that  the  steady  growth  of  frozen  soil  occurs  in  the 
region  bounded  by  curve  1  and  Ri .  Since  V0  is  an 
arbitrary  positive  number,  eq  69a  also  implies  that 
the  steady  growth  of  frozen  soil  occurs  every¬ 
where  in  Ry.  However,  the  steady  growth  of  ice- 
rich  frozen  soil  is  anticipated  to  occur  in  the  part  of 
R„  not  far  from  the  boundary  Ri . 

Suppose  that  a  point  (oq/io)  in  is  given;  then 
we  can  find  V0  that  satisfies  eq  69d.  At  the  end  of 
the  preceding  section  we  had  three  independent 
variables,  Oq,  oq  and  Vq.  Since  these  three  variables 
are  related  by  eq  69d,  we  now  have  only  two 
independent  variables,  Oq  and  cq.  In  other  words, 
we  have  found  that  the  condition  of  the  steady 
growth  of  ice-rich  frozen  soil  is  uniquely  deter¬ 
mined  by  two  independent  variables,  otg  and  cq. 
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under  a  given  hydraulic  condition  and  overbur¬ 
den  pressure.  We  have  shown  that  there  exists  a 
traveling  wave  solution  containing  two  indepen¬ 
dent  parameters,  Oq  and  at,  to  the  problem  of 
steadily  growing  ice-rich  frozen  soil. 

DISCUSSION 

One  of  die  outstanding  questions  among  re¬ 
searchers  of  frost  heave  has  been  the  relationship 
between  the  rate  of  frost  heave  r  and  the  speed  of 
a  frost  front  V0.  A  significant  amount  of  effort  has 
been  made  to  determine  empirically  this  relation¬ 
ship  under  the  hypothesis  that  r  is  uniquely  deter¬ 
mined  by  Vq.  The  empirical  relationships  between 
r  and  V0  reported  in  the  literature  sometime  dis¬ 
agree  with  (or  evencontradict)  each  other  (Takashi 
et  al.  1978),  and  there  appears  to  be  no  consensus 
among  researchers.  This  situation  casts  a  serious 
doubt  upon  the  validity  of  the  underlying  hypothe¬ 
sis  that  r  is  uniquely  determined  by  V0. 

We  will  show  that  r  is  not  uniquely  determined 
by  V0  if  Mj  and  the  assumptions  used  in  our 
analysis  are  valid.  Using  eq  64b,  we  will  reduce  eq 
34f  to 

dzr  =  IT1  (fai  ai  -  *o<*o) 

~P9oVi(d|"ld2Wo-vi).  (70a) 


Since  and  Oq  are  related  by  eq  69d,  we  can 
express  r  as  a  function  of  V0  and  either  Oooroq.We 
will  write  eq  70a  in  two  ways  as 


d2r  —  SjOq  +  S2 

(70b) 

=  S3  Ctj  +  S4 

(70c) 

where  Sj  is  given  as 

11 

if 

(70d) 

S2  *  (l  -  d  i 1  dz )  p3o  yowo 

(70e) 

S3  =  WC2I  *21  &21 

(700 

-{l  -d^dalviJpjoVo. 

St  a  -  [(d? dx  -  kob2) (two  -  Vi ) 

(70g) 

We  will  examine  the  value  of  S(  for  a  special  case  in 
which  o  vanishes  and  V0  is  much  less  than  /10.  The 
limiting  value  of  Ss  as  V0 approaches  zero  is  given 
as 


Si  =  b*KSi 

(71a) 

Sj  =  0 

(71b) 

S3- fr-JC^S 

(71c) 

& 

II 

0 

(71d) 

S=*2i(ko  +  Lb*K#. 

(71e) 

It  follows  from  eq  71a,  71b,  71c  and  71d  that  die  last 
term  on  the  right-hand  side  of  eq  70b  (or  70c)  is 
negligible  in  comparison  with  the  first  term  in  die 
right-hand  side  of  eq  7%  (or  70c)  when  V0  is  much 
less  than/10and  oq  2: 1.0 °C  cm-1  for  i  *  1,2.  It  is  easy 
to  find  that  Sx,  ^  and  S3  are  positive,  but  the  sign 
of  S4  is  not  certain. 

According  to  die  results  of  our  analysis,  we 
have  found  that  die  rate  of  frost  heave  r  is  not 
uniquely  determined  by  die  speed  cf  a  frost  front 
V0  alone  and  that  die  relationship  between  r  and 
V0  strongly  depends  on  04(1 =0,1).  To  examine  the 
validity  of  eq  70b  we  will  use  reported  experimen¬ 
tal  data.  Takashi  et  al.  (1978)  conducted  a  series  of 
frost  heave  tests  in  which  the  temperature  of  die 
unfrozen  part  Rq  was  kept  constant  at  0.2-03°C 
higher  than  die  freezing  point  of  a  sample  so  that 
the  speed  V0  of  a  frost  front  n0  was  kept  nearly 
constant.  Dividing  eq  70b  by  d2  Vq,  we  will  reduce 
eq70bto 

l  =  tV?  =  d?bK2iaotf 

+  {d?-d;')p3ow0.  (72) 

Takashi  etal.  (1978)  called  £  the  "frost  heave  ratio. " 

Analyzing  their  data,Takashiet  al.  (1978)  found 
empirically  that  |is  uniquely  determined  by  V0for 
a  given  applied  pressure  a.  A  typical  behavior  of  | 
vs.  V0  obtained  by  them  is  reproduced  in  Figure  4 
where  a  curve  is  drawn  that  approximately  repre¬ 
sents  their  data  points  taken  with  their  sample  2 
under  die  applied  pressure  a  =  304  kPa.  In  their 
tests  the  temperature  profile  in  the  sample  was  not 
measured  and  it  is  difficult  to  assess  the  variability 
of  Oq.  However,  if  Oq  is  kept  nearly  constant  and 
fC2i  mainly  depends  onao  VQ~\  then  their  empirical 
relationships  between  §  and  V0  are  consistent  with 
eq  72  as  we  will  show  below. 

Since  die  second  term  on  die  right  side  of  eq  72 
is  a  given  constant,  |  approaches  asymptotically 
this  constant  as  V0  becomes  infinite.  Thevalue  of  £ 
increases  with  the  decreasing  V0  until  £  becomes 
infinite  when  V0  vanishes  and  an  ke  layer  grows. 
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Figure  4.  Frost  heave  ratio  £  (%)  os.  V0  (cm  (T1) 
obtained  empirically  by  Takashi  et  al.  (1978). 

A 

The  curve  vs.  V0  in  Figure  4  is  converted  into  the 
curve  £  vs.  V(f 1  in  Figure  5.  It  is  clear  from  eq  72  that 
the  gradient  of  the  curve  in  Figure  5  is  proportional 
to  fC2i  if  eq  72  holds  true.  FromFigure5we  find  that 
the  gradient  of  the  curve  tends  to  decrease  with  the 
increasing  Vo"1;  namely,  K21  is  a  decreasing  func¬ 
tion  of  Vff1 . 

We  have  derived  eq  72  under  the  assumption 
that  the  speed  of  a  frost  front  V0  is  constant.  There¬ 
fore,  eq  72  is  not  anticipated  to  hold  true  for  the 
transient  freezing  in  which  V0  varies  with  time. 
However,  eq  72  may  approximately  hold  true  for 
the  transient  case  in  which  the  change  of  V0  with 
time  is  small.  Analyzing  the  data  on  transient 
freezing  tests  obtained  by  Akagawa  (1990),  Miyata 
and  Akagawa  (1991)  empirically  found  that  £  may 
be  uniquely  determined  by  a0^\  though  data  are 
limited.  Their  data  £  of  two  tests  (test  A  with  <3=60 
kPa  and  test  B  with  o = 110  kPa)  vs.  aoVJ1  are  pre¬ 
sented  in  Figure  6  where  a  curve  is  drawn  to  show 
the  trend  of  the  data  points. 

A  soil  specimen  of  85-cm  length  was  frozen 
from  the  bottom  up  with  constant  boundary  tem¬ 
peratures  (Akagawa  1990),  The  data  points  (£, 
do  Vo-1)  were  taken  during  the  time  period  from  4  to 
45  hours  after  the  start  of  the  test  The  speed  VQ 
decreased  and  <Xq  increased  monotonkally  with 
time;  hence,  the  value  OoVfT1  increased  monotoni- 
cally  with  time.  In  test  A,  for  instance,  V0  changed 
from  9.12  cm  d"1  at  4  hours  to  nearly  zero  at  45 
hours  while  Oq  changed  from  about  0.6°C  cm-1  at 
4  hours  to  about  1.5°C  cm-1  at  23  hours.  It  is  quite 
interesting  that  eq  72  may  hold  true  despite  such 


Figure  5.  Frost  beam  ratio  £  (%)  re.  Vq1  (cmr1  d) 
obtained  empirically  by  Takashi  et  al.  (1978). 
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Figure6.  Frost  heave  ratio  £  (%)  re.  ooVo1  (°Ccnr2d) 
obtained  empirically  by  Miyata  and  Akagawa  (1991). 

rates  of  change  in  V0  and  (Xq  as  described  above. 
From  Figure  6  we  find  a  trend  similar  to  that  of 
Figure  5:  the  gradient  of  foe  curve  £  vs.  OoV^1  tends 
to  decrease  with  the  increasing  OoVfr1,  namely,  K2i 
is  a  decreasing  function  of  oCoVfr1. 

We  have  studied  foe  steady  growth  of  ice-rich 
frozen  soil  by  using  Mj.  We  have  shown  that  there 
exists  a  traveling  wave  solution  to  foe  problem  of 
steadily  growing  ice-rich  frozen  soil  and  that  this 
solution  is  reduced  to  foe  solution  to  foe  problem 


of  a  steadily  growing  ice  layer  when  the  velocity 
no  vanishes.  We  have  also  shown  that  the  steady 
growth  condition  of  ice-rich  frozen  soil  under 
given  hydraulic  conditions  and  applied  pressures 
is  uniquely  determined  by  a  set  of  two  physical 
variables,  Oq  and  at.  We  will  present  the  results  of 
our  experimental  study  in  another  report. 
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APPENDIX  A:  EXACT  SOLUTION  OF  EQUATION  45 


When  *,(£)  is  given  by  eq  40e,  we  will  introduce  a 
new  independent  variable  X  defined  for  nt  S^S^as 

X=l+n(4-/io),  1+tlSiX^l.  (Al) 

Using  X,  we  will  reduce  eq  45  to 

=  (nX)'1  (PiT  -  txo  -  ifco1  AL) .  (A2) 

dX 

Multiplying  eq  A2  by  X-®1  we  will  write  eq  A2  as 

J_  TX-(Pi  A)  =  _ (^xjT1  X'(p,/4) (a^AL) . 
dX 

(A3) 

Suppose  that  a  solution  7T(X)  of  eq  A2  exists.  Sincejhe 
right  side  of  eq  A3  is  negative,  the  function  TX-*1  **  is 
a  decreasing  function  of  X.  Therefore,  7\X)  must  be 
negative  forX>  1  because  T=0atX=  1.  Wehave  found 


that  a  solution  1\X)  is  negative  for  X>  1  if  it  exists. 
Integrating  eq  A2,  we  obtain 

IT 

vifaT-ao-tf  A iX'dT.  (A4) 

o 

Since  the  integrand  of  eq  A4  is  continuous,  a  solution 
7\X)  exists.  It  is  easy  to  see  that  the  right  side  of  eq  A2 
satisfies  a  Lipschitz  condition  with  respect  to  Tbecause 
the  function  A(T)  possesses  a  continuous  first  deriva¬ 
tive  as  assumed.  Based  on  an  elementary  theorem  of 
ordinary  differential  equations  (Sansone  and  Conti 
1964),  we  may  conclude  that  eq  A2  has  a  unique 
negative  solution  1\X)  forX  >  1. 

It  follows  from  eq  A2  that  the  unique  solution  of 
eq  A2  is  decreasing  with  x  (or  4)  because  the  right 
side  of  eq  A2  is  negative.  We  have  shown  that  eq  45 
has  a  unique  and  decreasing  solution  for  nt  £  £  ^  Hq. 
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APPENDIX  B:  APPROXIMATE  SOLUTION  OF  EQUATION  45 


As  we  have  shown  above  (App.  A),  eq  45  has  a 
unique  and  decreasing  solution.  This  implies  that  X  (or 
4)  and  T  are  one-to-one.  Treating  A  as  a  function  of  X, 
we  will  write  eq  A2  as 

T(X)  =  mffi1  (l  -  XM'1)  +  non-1  YX*m'X  (Bl) 
where  Xq,  X|,  and  Kare  defined  as 


*0  =  kq~*Lp3o  Vo 

(B2) 

xi  =  oco  +  no  wo 

(B3) 

fX 

Y(X)  =  j  vX'^^iX. 

(B4) 

We  will  rewrite  eq  A2  as 

ir.-ttxr'mYi 

(B5) 

where  Yx  is  defined  as 

Yi  =  1  -  Xi1  (Pi  T  +  no  v). 

(B6) 

Using  eq  B5,  we  will  reduce  eq  B4  to 

Y  =  qxj1  v 

(B7) 

v  =-|  vX'Pm\tdT. 

(B8) 

Using  eq  B7,  we  will  reduce  eq  Bl  to: 

T(4)  =  m pi1  (l  -  XM_1 )  +  icoxI1vXm‘1  .  (B9) 


Using  eq  B5  and  B9,  we  obtain: 

7*(4)X  =  -  X^Ui  -  xoni1  ft v)  +  xov.  (BIO) 


We  will  seek  approximate  equations  for  eq  B9  and 
B 10  when  the  following  condition  holds  true: 

1)8  <  1  and  Pt5<l.  (BU) 

We  found  that  the  condition,  eq  BU.  holds  true  when 
the  steady  growth  of  an  ice  layer  occurs  (Takeda  and 
Nakano  1990).  When  eq  B 1 1  holds  true,  we  may  use  an 
approximation  (Nakano  1990)  given  as: 

X***  =l  +  pi(5-n0).  (B12) 

Using  eq  B12,  we  will  reduce  eq  B9  and  BIO  to 

7(4)  =  Roxi1  v  - (*i -  RoxJ1  v  pi)(§  -  «o)(B13) 

T(f,)X  =-m  +  7Co  (v  +  xj1  Pi  v) 

-  P»(*i  -  itoxi*  Pi v)  (4  -  no).  (B14) 

For  the  calculation  of  v,  X  *  is  approximated 
by  using  eqB  12  as 

X^1”'1  =  [1  +  Pi  (£  -  no)]"1.  (B15) 

When  V0  is  very  small,  from  eq  B13  T(4)  *s  given  as 
7X4)  =  -«i(4-«o)-  (B16) 

Substituting  eq  B16  into  eq  B15  we  obtain 

X-*n-1  =  (!_  plJl-1T)"1.  (B17) 

Now  v  is  approximated  as 

,T 

v  =  -  j  vIYi  (1  -  Pi  Ki'TjfdT.  (B18) 
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APPENDIX  C:  COMPUTATION  OF  I, 


Using  eq  B5,  we  will  reduce  eq  56b  to 

-niKio h  =  f  X(Kw /Ki)YildT.  (Cl) 

Jo 

The  term  P1TineqB6  describes  the  effect  of  sensi¬ 
tive  heat  that  is  much  less  than  xt.  Also  the  term 
xlSco  v  in  eq  B6  is  less  than  one.  Therefore,  the  term 
KlHPi T  +  xov)  is  generally  less  than  one  and  we 
may  approximate  Yj1  as 

Yt^l  +  xftPiT  +  Kov).  (C2) 

Using  eq  55b,  55c  and  C2,  we  will  reduce  eq  Cl  to: 


-  «i  Xio /i  -  a>0  no  w^Ti  ^2i 

=  |  ♦iT,(l+*I1pi7VT 

=  ♦nTi(l  -e)  (C3) 

where 

♦i  =  d9i/dT 

e-M^j-Ti  +  fuj  ♦idr|  (C4) 
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